Theoretical study of the chlorine adsorption on the Ag(lll) surface 
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We study the adsorption of chlorine on the Ag(lll) surface with full potential gradient corrected 
density functional calculations. When considering a y3 x v3 R30° pattern, we find that the fee 
hollow is the most favorable adsorption site. We obtain an Ag-Cl bond length of 2.62 A which is 
intermediate between two controversial experimental results. We discuss the differences of core level 
energies and densities of states for the different adsorption sites. We find the Cl-Ag interaction to 
be more consistent with an ionic, rather than covalent, picture of the bonding. In addition, we 
compute energetics and related properties of Ag bulk and the clean Ag surface comparing the local 
density approximation and the generalized gradient approximation. 



o 
o 

(N 



o 



in 



X 



I. INTRODUCTION 



The study of halogenide adsorption on metal surfaces 
has become a very interesting field in surface science. 

S Important questions such as adsorption geometries and 
energetics have been addressed in experimental stud- 
'-^J \ ies. The adsorption of chlorine on the Ag(lll) sur- 
£5 ■ face has been studied by various groups and with dif- 
ferent techniques including low energy electron diffrac- 
i tion (LEED)ETEi, surface extended X-ray absorption fine 

structure (jSEXAFS)BO, Auger spectroscopytrEl, thermal 
7— I , desorptioia rU, measurements of adsorption isobars and 
isothermsB work function measurements!! and scanning 
tunneling microscopyu'lj. 



However, there is still an ongoing discussion on the 
correct interpretation of the experiments. Various expla- 
nations of observed adsorption patterns have been pro- 
posed, and the bond length has, been measured twice with 
significantly different resultsStj. 

We recently presented a theoreticaL-study of the chlo- 
rine adsorption on the Cu(lll) surfaced and we are there- 
fore confident that simulation can shed some light on the 
difficult question of chlorine adsorption. As a test of our 
computational method, we performed additionally cal- 
culations on Ag bulk and the clean surfaces. These are 
presented in section ||. The question of chlorine adsorp- 



tion is discussed in section [H 



II. AG BULK AND CLEAN SURFACES 
A. Computational Parameters 



We used a local basis set formalism where the ba- 
sis functions are Gaussian type orbitals centered at the 
atoms as implemented in the code CRYSTALEJ. In or- 
der to examine the effect of various treatments of electron 
exchange and correlation, we have used the local density 
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approximation (LDA) with Dirac-Slater exchang and 
the Perdew-Zunger correlation functional^, the gradient 
corrected exchange aneLcorrelation functional of Perdew 
and Wang (GGA-PW)E3 and the gradient corrected func- 
tional of Perdew, Burke and Ernzerhof (GGA-PBE)EJ. 

For silver, we used ajjelativistic pseudopotential with 
corresponding basis setti For the calculation on the free 
atom, we used the original basis set. For the calcula- 
tion on the solid and surfaces, we kept the inner expo- 
nents ([2slpld]) fixed and re-optimized the outermost ex- 
ponents within the GGA-PW approximation. Values of 
0.63 4r (with the Bohr radius ao) and 0.11 -\ for the s- 

a Q a 

exponents, 0.23 \ for the p-exponent and 0.18 \ for the 

a Q a Q 

d-exponent were obtained (the most diffuse s-exponent 
could not be optimized and was kept fixed at a value of 
0.11 -tt to avoid numerical instabilities). This resulted 

as a whole in a [4s3p2d] basis set. The CI basis set is 
[5s4pld]_as used in our recent study of CI adsorption on 
Cu(lll)H. 

The density functional potential was fitted with an 
auxiliary basis set, which consisted, for both Ag and CI, 
of 12 s and p functions with exponents taken to be a ge- 
ometrical sequence from 0.1 ^ to 2000 For CI 5 d 

functions with exponents in the range 0.8 A- to 100 \ 

were employed, for Ag 8 d and / functions from 0.1 

to 100 . 

We used a Pack-Monkhorst net of shrinking factor 16 
for Brillouin zone sampling, which resulted in 145 irre- 
ducible fc-points for the bulk, 45 for the clean (100) sur- 
face, 81 for the clean (110) surface, 30 for the clean (111) 
surface, and 30 or 73 for the system Cl/Ag(lll), depend- 
ing on the adsorption site (30 for hep, 73 for fee, bridge 
and atop site — an explanation of the sites is given in sec- 
tion |n|; see a l so figure [l 13) . A smearing temperature of 
0.03 E h was applied to the Fermi function (£^=27.2114 



eV) . This approach has been 
for lithium bulk and surface.' 
of CI on the Cu(lll) surface! 



3wn to give stable results 
LZl as well as for the study 



B. Properties of bulk Ag and clean Ag surfaces 

Results for cohesive properties of bulk Ag are in ex- 
cellent agreement with those published in the literature 
(table Q). The tendency of the LDA to overbind results in 
a short lattice constant, too high bulk modulus and too 
high binding energy. GGA improves the binding energy 
and bulk modulus; the deviation of the lattice constant 
is similar to LDA (with opposite sign, however). 

Surface formation energies (in Eh per surface atom) 
of the clean, unrelaxed surfaces are displayed in table |J. 
The surfaces were modeled with a finite number of Ag 
layers (typically from 1 to 7); this model was not period- 
ically repeated in space. The surface energy was found 
to be stable from three layers onwards. The bulk energy 
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was computed with the same computational parameters 
(basis set, shrinking factors, smearing temperature). 

LDA binding energies are approximatively 50% higher 
than those obtained with gradient corrected fusctionals, 
which is higher than found in a previous studyfci GGA- 
PW and GGA-PBE turned out to give practically iden- 
tical results in all cases which are in excellent agreement 
with the literature!!!]. 



III. CL ON THE AG(lll) SURFACE 

The adsorption of CI on the Ag(lll) surface has been 
investigated by several groups. Still, there is no agree- 
ment on adsorption patterns, coverages, preferred sites 
and bond lengths. Jn the first study with LEED and 
Auger spectroscopy^, a \/~3 x V3 R30° pattern at low 
chlorine coverage and a (3 x 3) pattern at higher coverage 
were observed. In a second study, a v^x V3 R30° pattern 
at a coverage of 1/3 ancLa (10 x 10) pattern at a cov- 
erage of 0.49 were foundtl. Within an Arrhenius model, 
the binding energy was determined to be 0.080 E/UCh ■ 
A similar energy was found in another experiment!], but 
the LEED pattern was interpreted as an epitaxial AgCl 
layer. A -\/3 x \/3 R30° pattern at a coverage close to one 
was found in ref. |^. Scanning tunneling microscopy ex- 
periments have been interpreted as a parallel double row 
structure, where CI atoms occupy bridge and threefold 
hollow sitesa, or as (17 x 17) overlayer at saturationQ. 

A structural determination with SEXAFS gave a Ag- 
Cl bond length of 2.70 ± 0.02 A at a nominal coverage 
of 1/3 of a monolayer and 2.71) ± 0.01 A at a nominal 
coverage of 2/3 of a monolayers. CI was found to adsorb 
in threefold hollow sites in a y/3 x \/3 R30° pattern, for 
the higher coverage a vacancy honeycomb structure was 
proposed. This determination has been questioned in 
a recent article where a bond length of 2.48±0.04 A has 
been deduced from SEXAFS measurements for a nominal 
coverage of 1/3 of a monolayers. In this study an order- 
disorder transition was found at 195 K, below which a 
V3x \/3 R30° pattern at a coverage of 1/3 of a monolayer 
was identified in LEED measurements; this pattern was 
found to vanish above the transition temperature. In 
addition, a (13 x 13) overlayer was identified at a coverage 
of 0.37. The bond length of 2.48 A was found to fit better 
with other bond lengths of chlorine adsorbates on 3c? and 
4d fee metals. Therefore, it was suggested that previous 
experiments should be reinterpreted. A variation of the 
bond length with coverage was neither found for GLon 
Ag(lll)l3 nor for the similar system of I on Rh(lll)tj. 

These experimental observations therefore pose an in- 
teresting problem for simulation studies. In our calcula- 
tions, we modeled the adsorption of CI on Ag(lll) as a 
slab periodic in two dimensions but with a finite thick- 
ness of three or four layers of silver. This was found to be 
sufficiently thick to obtain stable results for the surface 
energy of the clean Ag(lll) surface. The slab was cov- 
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ered by a layer of chlorine corresponding to a coverage 
of one third of a monolayer, in a y/3 x V3 R30° pattern. 
CI was adsorbed in four different sites: the atop site on 
top of a silver atom in the top layer, the bridge site in 
the middle of two Ag in the top layer, and hep (and fee) 
hollows where the chlorine sits in threefold hollow sites 
with a silver atom in the second (third) layer under the 
chlorine atom (see figure [l]) . In the following, all the cal- 
culations were performed at the GGA-PW level, if not 
explicitly stated otherwise. 

The results of these calculations are summarized in ta- 
ble [II. First, it turns out that the threefold hollow sites 
are lowest in energy, with the bridge site being about 
0.003 Eh higher and the atop site 0.017 Eh higher in en- 
ergy. This picture is quite similar to that of chlorine ad- 
sorption on the Cu(lll) surfaces. As the threefold sites 
are nearly degenerate in energy, we considered an im- 
proved model with four silver layers to investigate the de- 
pendence of the results with respect to the number of lay- 
ers. It turns out that the fee hollow stays lower in energy, 
although the energy difference is reduced from 0.001 Eh 
(three layers) to only 0.0004 Eh (four layers). Therefore, 
a fully reliable prediction of the adsorption site can not be 
made as more layers might still change this picture. The 
total adsorption energy is 0.111 Eh per chlorine atom. 
We computed a dissociation energy Dq — 0.099 Eh, vi- 
brational constant lo p = 520— and internuclear distance 

L cm 

r e = 2.05 A for the chlorine dimer which is in reasonable 
agreement with experiment (Dq=0.091 Eh, uj e — 560—, 
r e =1.99 A)§. Therefore, our binding energy corresponds 
to 0.123 Eh per CI2 when computed with respect to the 
CI2 molecule, using our computed binding energy for the 
chlorine dimer. This is higher_than that deduced from 
thermal desorption experimentally. 080 E^V-and measure- 
ments of adsorption isobars and isothermsQ (~ 0.088 Eh 
for the heat of adsorption of CI2 ). However, the first 
experimental value depends on an empirical prefactor as- 
sumed in the Arrhenius model and was questioned in ref. 
|l|; the authors stated that it should be higher as Dq of 
the CI2 molecule because desorption as atomic CI was 
observed. 

Another important aspect is the adsorption geometry. 
We allowed the top silver layer to relax vertically and a 
slight inwards relaxation by about 0.04 A was found for 
all adsorption sites. The Ag-Cl bond length was deter- 
mined to be 2.62 A for the threefold hollow sites, 2.54 
A for the bridge site and 2.38 A for the atop site. Tests 
on the 4-layer slab demonstrated that these bond lengths 
are stable with respect to the number of layers. Other 
tests such as varying the outermost diffuse exponents in 
the basis set or applying a smaller smearing tempera- 
ture were performed, and it turned out that the bond 
length was also stable with respect to these parame- 
ters; the variation was of the order ± 0.02 A. There- 
fore, we can conclude that the computed bond length is 
not highly sensitive to computational details. The re- 
maining approximation in our simulations is the choice 
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of the exchange-correlation functional. The effects of this 
choice were examined by recomputing the surface prop- 
erties within the LDA (table IV). LDA gives the same 
preference of adsorption sites and a similar energy split- 
ting. The adsorption energies are 20 - 25 mEh higher 
than those computed using the GGA-PW, and the bond 
lengths are shorter by 0.05 - 0.07 A. The bond length 
computed with our model (2.62 A in GGA-PW, 2.55 A 
in LDA) is intermediate between two controversial ex- 
periments (2.48 AD and 2.70 AO); therefore we cannot 
support either experiment over the other. 

Combining our data for CI on Ag(lll) and CI on 
Cu(lll), we can comment on the effective radius of the CI 
adsorbate, under the assumption that both the chlorine 
adsorbate and the copper and silver surfaces are made of 
hard spheres. For the fee site and CI on Cu, a bond length 
of 2.40 A was obtained. If we subtract the atomic radius 
of copper of §^7§A, then we obtain a radius for CI of 1.12 

A. For CI on Ag, we obtain a radius of 2.62A-^2A= 

1.17 A. This radius for CI is intermediate-between the 
atomic (0.99 A) and ionic (CI - ; 1.81 A))EI. The bond 
lengths predicted here and previously computed for the 
Cu(lll)-(-\/3 x y3) R30°-C1 give a consistent picture of 
the radius of the adsorbed Cl~ ion. 

The Mulliken population analysis gave a slight charge 
transfer of 0.2 e to 0.25 |e| to the chlorine atom, increas- 
ing in the order fee ~ hep < bridge < atop ; of course, 
one has to keep in mind that the error associated with 
this number is larger in the case of metals with diffuse 
basis functions because the diffuse CI basis functions can 
be used to partially describe Ag states and vice versa. 
Still, it is very interesting to compare the charge transfer 
with the position of the CI 3s eigenvalue (table |y|) . If we 
assume that intra-atomic interactions dominate the core 
eigenvalues of the CI atom, then the position of the CI 3s 
eigenvalue provides a measure of the atomic charge: the 
more negatively charged the CI ion, the lower the binding 
energy of the core eigenvalue. As this is related to the 
observable core level shift, it can provide an experimen- 
tal mechanism for distinguishing adsorption sites. We 
observe this correlation in our results for different ad- 
sorption sites (table |v|) , but we find that this correlation 
is not retained when comparing CI adsorption on Cu to 
that on Ag. The 3s eigenvalue is at lower binding en- 
ergy for CI on the Ag surface despite the CI charge being 
lower in magnitude than that of CI on the Cu surface. 
This could be due to errors in the Mulliken charge esti- 
mates which are more reliable within a consistent set of 
basis sets. 

The position of the CI levels is visible in the Mulliken 
projected densities of states (figures HI and |). First, we 
projected the density of states on the chlorine orbitals 
(figure 0)' The 3s level moves as described to higher 
energies (corresponding to lower binding energies) with 
higher Mulliken orbital populations. For all adsorption 
sites, the valence band contribution consists of a broad 
background from ~ — 0.2.E/J to ~ — Q.08Eh which con- 
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tains contributions from CI 3p x ,p y and 3p z . The peak 
at ~ —0.1 Eh is due to 3p x and 3p y contributions in all 
cases; there is no distinct peak from 3p z contributions. 
The 3p x peak and 3p y peak are virtually identical in the 
case of fee, hep or atop adsorption sites; for the bridge 
site the p y peak carries more weight than the p x peak 
(figure ||). These findings are consistent with orbital pro- 
jected Mulliken charges (table |Vl|) which are highest for 
p x and p y in the case of fee, hep and atop site. For the 
bridge site, the p y population is higher than p x and again 
higher than p z — note that, in our choice of geometry, 
for the bridge site, the CI 3p x orbital overlaps more with 
the silver atoms in the top layer than the CI 3p y orbital 
(see figure |l|). We can therefore conclude that in all sites 
the charge transferred to the CI atom occupies those or- 
bitals which have lowest overlap with the silver atoms. 
This is consistent with an ionic rather than covalent pic- 
ture of the bonding. We note that these orbitals may 
be separately observed by exploiting the polarization de- 
pendence of angle resolved photoemissionEj. When com- 
paring the chlorine adsorption on silver with the chlorine 
adsorption on copper (figure ^), we find that the band- 
width of the 3p states is similar; as mentioned earlier the 
3s peak is at slightly higher energy for CI on Ag than for 
CI on Cu, although the charge is a bit smaller in magni- 
tude (-0.207 |e| for copper vs -0.198 |e| for silver). The 
peak structure is different from CI on Ag(lll) as we find 
one additional peak for CI 3p z ; the two other peaks carry 
equal weight from CI 3p x and CI 3p y . As in the case of 
silver, the CI 3p x and 3p y orbitals carry slightly more 
charge than CI 3p z . 

Finally, in figures || and |], we display the total charge 
density of CI adsorbed in the fee hollow site and the 
charge density difference. The plane was chosen in such 
a way that the center of the chlorine atom, the center of 
one of the nearest neighbor silver atoms (at a distance 
of 2.62 A) and the centers of two more silver atoms (at 
a distance of 3.91 A) lie in the same plane. The charge 
density difference has been obtained by subtracting the 
charge densities of a clear Ag slab (with the geometry of 
the system Ag(lll)-(%/3 x R30°-C1 ) and a layer of 
chlorine atoms from the charge density of the system af- 
ter the chlorine adsorption. We find no evidence of bond 
formation in the charge density plots so that these plots 
support the idea of an ionic binding mechanism. A fur- 
ther evidence of an ionic bonding is the very small overlap 
population of ~ 0.01 |e| between CI and the nearest silver 
atom. 



IV. SUMMARY 

We have studied the adsorption of chlorine on the 
Ag(lll) surface using density functional theory calcula- 
tions in which we have paid close attention to numerical 
accuracy. The fee hollow site was found to be the pre- 
ferred adsorption site, slightly lower in energy than the 
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hep site. The adsorption energy is computed to be 0.111 
Eh per atom with respect to free CI atoms or 0.123 Eh 
per CI2 molecule. 

The computed Ag-Cl distance of 2.62 A at the GGA- 
PW level (or 2.55 A at the LDA level) arc intermediate 
between two controversial values deduced from experi- 
ment. The bond length is consistent with that computed 
for CI on Cu which would indicate an effective radius 
of the CI adsorbate of <~ 1.15 A. For the different ad- 
sorption sites of CI on Ag(lll), we find that the charge 
on chlorine is consistent with the relative position of the 
CI 3s level, i.e. the higher the charge on CI, the lower 
the binding energy of the 3s level. This might be valu- 
able in the interpretation of core level shift observations. 
From projected densities of states, we conclude that the 
chlorine charge is larger in those orbitals which are most 
separated from silver which means that the character of 
the bond tends towards ionic rather than covalcnt. In 
our choice of geometry, the chlorine 3p x and 3p y orbitals 
have a slightly higher occupation than the chlorine 3p z 
orbital for the threefold hollow sites which might also be 
observable in photoemission experiments. 
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TABLE I. The ground state properties of bulk Ag. 



LDA 

GGA-PW, GGA-PBE 



LDAi 
LDA@p 

GGA-PW^ 
GGA-PWEi 
experiment 



a [A] 
3.98 
4.10 

3.95 ... 4.01 

4.00 
4.12 ... 4.14 

4.03& 



Ecoh [Eh] 

0.129 
0.088 



0.109 



B [GPa] 
165 
113 

137 ... 163 

139 
93 ... 98 
84|l 

102 







TABLE II. The surface energy ( — e Efl ? — ) of the low index silver surfaces. 

0i/ ^ surracc atom ' 



surface 


LDA 


GGA-PW, GGA-PBE 


Ref. [[! (GGA-PBE) 


(100) 


0.033 


0.023 


0.024 


(110) 


0.050 


0.033 


0.035 


(111) 


0.026 


0.017 


0.020 
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TABLE III. Adsorption of CI on the Ag(lll) surface, GGA-PW results. <5i_ 2 is the change in interlayer spacing between 
first and second silver layers relative to the bulk value, dci-Ag top layer is the interlayer distance between the CI layer and 
the top Ag layer, dci-Ag nn is the distance between CI and nearest neighbor Ag. The adsorption energy is the difference 
Eci at Ag(in) — ^Agfin) — -E-C1- n is the number of silver layers in the slab. 



n 

fee site 


rfd-Ag top layer, in A 


Si-2, in A 


dci-Ag nn, in A 


Eadsorption in Eh, per CI atom 


3 


2.01 


-0.043 


2.62 


-0.11164 


4 


2.01 


-0.039 


2.62 


-0.11173 


hep site 
3 


2.01 


-0.042 


2.62 


-0.11068 


4 


2.02 


-0.036 


2.62 


-0.11129 


bridge site 
3 


2.08 


-0.040 


2.54 


-0.10878 


atop site 
3 


2.38 


-0.034 


2.38 


-0.09460 
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TABLE IV. Adsorption of CI on the Ag(lll) surface, LDA results. The notation is the same as in table LT 



site dci-Ag top layer, in A S1-2, in A d C i-A g nn, in A Eadaorption in E h , per CI atom 

fee 1.96 -0.052 2.55 -0.13653 

hep 1.96 -0.046 2.55 -0.13422 

bridge 2.03 -0.052 2.47 -0.13297 

atop 2.33 -0.047 2.33 -0.11487 
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TABLE V. Charge and position of 3s eigenvalue for CI on different adsorption sites. 



site 


charge, in e 


3s level relative to i?Fermi, in E h 




CI on Ag: 




fee 


-0.198 


-0.563 


hep 


-0.204 


-0.562 


bridge 


-0.218 


-0.555 


atop 


-0.252 


-0.532 




CI on Cu: 




fee 


-0.207 


-0.578 
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TABLE VI. Orbital projected CI charge, in e|, on different adsorption sites. 



site 


Px 


Py 


Pz 




CI on Ag: 






r 

ICC 


O 7H A 

-6. /94 


O 7H A 

-6. /94 


-3.62y 


hep 


-3.795 


-3.795 


-3.636 


bridge 


-3.766 


-3.840 


-3.628 


atop 


-3.822 


-3.822 


-3.605 




CI on Cu: 






fee 


-3.802 


-3.802 


-3.598 
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FIG. 1. The structures considered for CI, adsorbed on the Ag(lll) surface at one third coverage, in a 73 x ^3 R30° unit 
cell. When CI is adsorbed in an fee hollow, it sits above a Ag atom in the third layer (upper left) while in the hep hollow it is 
above an atom in the second layer (upper right). The atop position is CI adsorbed vertically above a surface atom (lower left). 
In the bridge position it is vertically above the middle of two surface atoms (lower right). The 2>direction is horizontal in the 
plane, the {/-direction is vertical, and the z direction points out of the plane of the page. 




atop bridge 
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FIG. 2. Projected density of states for CI on Ag(lll) and different adsorption sites. The part of the density of states which 
is due to the CI basis functions, is shown in this projection. The Fermi energy is fixed at Eh- 
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FIG. 3. Projected density of states for CI on Ag(lll), bridge site. The density of states is shown, projected on the CI basis 
functions (px,p y and p z ). The Fermi energy is fixed at -Eh. 
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FIG. 4. Projected density of states for CI on Cu(lll) and CI on Ag(lll), fee hollow. The part of the density of states which 
is due to the CI basis functions, is shown in this projection. The Fermi energy is fixed at Eh- 
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FIG. 5. Total charge density for CI on Ag(lll), fee hollow, from 0.01 to 0.09 a.u. in steps of 0.01. 
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FIG. 6. Difference in charge density (i.e. charge CI on Ag(lll) - charge of clear Ag(lll) - charge of a layer of CI atoms) 
for CI on Ag(lll), fee hollow, from -0.007 to 0.009 in steps of 0.001 a.u. Excess negative charge is displayed with long dashed 
lines, excess positive charge with solid lines, zero difference with short dashed lines 
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